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Abstract. It is shown that if the self-gravitating shearing sheet, a model of a patch of a galactic disk, is embedded 
in a live dark halo, this has a strong effect on the dynamics of density waves in the sheet. I describe how the 
density waves and the halo interact via halo particles either on orbits in resonance with the wave or on non- 
resonant orbits. Contrary to expectation the presence of the halo leads to a very considerable enhancement of 
the amplitudes of the density waves in the shearing sheet. This effect appears to be the equivalent of the recently 
reported enhanced growth of bars in numerically simulated stellar disks embedded in live dark halos. Finally I 
discuss the transfer of linear momentum from a density wave in the sheet to the halo and show that it is mediated 
only by halo particles on resonant orbits. 
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1. Introduction 

The shearing sheet (Goldreich & Lynden-Bcll 1965, Julian 
& Toomre 1966) model has been developed as a tool to 
study the dynamics of galactic disks and is particularly 
well suited to describe theoretically the dynamical mech- 
anisms responsible for the formation of spiral arms. For 
the sake of simplicity the model describes only the dynam- 
ics of a patch of a galactic disk. It is assumed to be in- 
finitesimally thin and its radial size is assumed to be much 
smaller than the disk. Polar coordinates can be therefore 
rectified to pseudo-Cartesian coordinates and the velocity 
field of the differential rotation of the disk can be approx- 
imated by a linear shear flow. These simplifications allow 
an analytical treatment of the problem, which helps to 
clarify the underlying physical processes operating in the 
disk. 

In the present paper of this series I discuss the dynam- 
ical effects if the shearing sheet is immersed in a live dark 
halo. Dark halos are usually thought to stabilize galac- 
tic disks against non-axisymmetric instabilities. This was 
first proposed by Ostriker & Peebles (1973) on the basis 
of - low-resolution - numerical simulations. Their phys- 
ical argument was that the presence of a dark halo re- 
duces the destabilizing self-gravity of the disks. Doubts 
about an entirely passive role of dark halos were raised 
by Toomre (1977), but he (Toomre 1981) also pointed out 
that a dense core of a dark halo may cut the feed-back loop 



of the corotation amplifier of bars or spiral density waves 
and suppress thus their growth. Recent high-resolution nu- 
merical simulations by Debattista & Sellwood (2000) and 
Athanassoula (2002, 2003) have shown that quite the re- 
verse, a destabilization of disks immersed in dark halos, 
might be actually true. Athanassoula (2002) demonstrated 
clearly that much stronger bars grow in the simulations if 
the disk is embedded in a live dark halo instead of a static 
halo potential. This is attributed to angular momentum 
transfer from the bar to the halo via halo particles on res- 
onant orbits. Angular momentum exchange between disk 
and halo has been addressed since the pioneering work of 
Weinberg (1985) in many studies theoretically or by nu- 
merical simulations and I refer to Athanassoula (2003) for 
an overview of the literature. Toomre (1981) has shown 
how the bar instability can be understood as an interfer- 
ence of spiral density waves in a resonance cavity between 
the corotation amplifier and an inner reflector (cf. also 
Fuchs 2004). Thus it is to be expected that a live dark 
halo will be also responsive to spiral density waves and 
develop wakes, which I investigate here using the shearing 
sheet model. 

2. Boltzmann equations 

The evolution of the distribution function of the disk stars 
in phase space is described by the linearized 4— dimensional 
Boltzmann equation 



Correspondence to: fuchs@ari.uni-heidelberg.de 



dfdl u dfdl 
dt dx 



dy 



2 



B. Fuchs: Density waves in the shearing sheet 



<9$do + $ho dfdi 9$do + $ho df d i 



dx du 
3$ d i + $ hl df d0 



dy dv 

d$ dl + $ hl a/ d0 



9x 



9u 



0, 



(1) 



where (x, y) denote spatial coordinates with y pointing in 
the direction of galactic rotation and (u, v) are the cor- 
responding velocity components, respectively. Equation 
(1) has been derived from the general 6 dimensional 
Boltzmann equation assuming delta function like depen- 
dencies of the distribution function on the vertical z coor- 
dinate and the vertical w velocity component, respectively, 
and integrating the Boltzmann equation with respect to 
them. A perturbation Ansatz of the form 

/d = /do + /di , $d = $do + $di , $h = $ho + $hi (2) 

is chosen for the distribution function of the disk stars 
and the gravitational potentials of the disk and the halo, 
respectively, and the Boltzmann equation (1) has been 
linearized accordingly. 

Similarly the linearized Boltzmann equation for the 
halo particles can be written as 
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(3) 



The choice of the dark halo model is lead by the following 
considerations. One of the deeper reasons for the success of 
the infinite shearing sheet model to describe spiral density 
waves realistically is the rapid convergence of the Poisson 
integral in self-gravitating disks (Julian & Toomre 1966). 
Consider, for example, the potential of a sinusoidal density 
perturbation 
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where G denotes the constant of gravitation. At x = 

$ = -4GE 10 sin(fc 2 ;) lim = 

27rGSio sin (kx) 
k ' 



,(4) 



(5) 



The sine integral in equation (5) converges so rapidly that 
it reaches at kx^ = § already 87% of its asymptotic value. 
Thus the "effective range" of gravity is about only a quar- 
ter of a wave length. The shearing sheet models effectively 
patches of galactic disks of such size. The wave lengths of 
density waves are of the order of the critical wave length 
(Julian & Toomre 1966) 
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where k denotes the epicyclic frequency of the stellar or- 
bits and Ed is the surface density of the sheet. In the 
solar neighbourhood in the Milky Way, for instance, the 
critical wave length is A cr j t = 5 kpc. Thus it is reason- 
able to neglect over such length scales, like in the shearing 
sheet model, the curvature of the mean circular orbits of 
the stars around the galactic center or the gradient of the 
surface density. The curvature of the stellar orbits due 
to the epicyclic motions of the stars, on the other hand, 
cannot be neglected and is indeed not neglected in the 
shearing sheet model. The radial size of an epicycle is ap- 
proximately given by cr u / 'k, where er u denotes the radial 
velocity dispersion of the stars, and the ratio of epicycle 
size and critical wave length is given by 



= 0.085Q 



(7) 



in terms of the Toomre stability parameter Q (Toomre 
1964) which is typically of the order of 1 to 2. Concurrent 
to these approximations I assume a dark halo which is 
homogeneous in its unperturbed state. Accordingly the 
curvature of the unperturbed orbits of the halo parti- 
cles is neglected on the scales considered here and the 
particles are assumed to be on straight-line orbits. The 
equations of motion of the halo particles are the charac- 
teristics of the Boltzmann equation. In a homogeneous 
halo x = V($do + ^ho) = 0, and in accordance with 
this assumption I neglect the force terms V<I>do and V$ho 
in the Boltzmann equation (3). This simplifies its solu- 
tion considerably. The disadvantage of such a model is 
that there arc no higher-order resonances of the orbits of 
the halo particles with the density waves as described by 
Weinberg (1985) or observed in the high-resolution simu- 
lations by Athanassoula (2002, 2003). However their effect 
was shown to be much less important than the main reso- 
nances of the particles with the density waves, which are 
properly described in the present model. 



3. Halo dynamics 

The Boltzmann equation (3) can be viewed as a linear 
partial differential equation for the perturbation of the 
distribution function of the halo particles, /hi, with inho- 
mogeneities depending on the perturbations of the grav- 
itational potentials of the disk and the halo, V$di and 
V$hi, respectively. Thus the equation can be solved for 
the disk and halo inhomogencitics separately and the so- 
lutions combined afterwards by superposition. 

3.1. Jeans instability of the halo 

I consider first that part of the Boltzmann equation (3) 



dfhi u dfh\ v dfhi tt , d/hi 
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Equation (8) is Fourier transformed with respect to time 
and spatial coordinates and I consider in the following 
Fourier terms 



/hk.u,, *hk,^ x e 



;(wt+(k,x)) 



(9) 



with frequency u) and wave vector k. If orthogonal spatial 
coordinates £, rj, and £ with £ parallel to the wave vector 
k are introduced, equation (8) takes the form 



hO 



= 0. 



(10) 



where k = |k| and v denotes the velocity component par- 
allel to the £-axis. Equation (10) can be immediately in- 
tegrated with respect to the velocity components perpen- 
dicular to k. Adopting for the background distribution 
function a Gaussian distribution, /ho = 
the solution of equation (10) is given by 
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The corresponding Fourier coefhcicnt of the density per- 
turbation of the halo particles is 



phk = -'I'hk 



Pb 



dv- 



The perturbations of the halo are supposed to be self- 
gravitating and the gravitational potential has to solve 
the Poisson equation 
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Combining equations (12) and (13) leads to the dispersion 
relation 
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which is well known from plasma physics (cf. Kegel 1998). 
Indeed, the dispersion relation (14) describes simply the 
Jeans collapse of the halo. From the imaginary part of 
equation (14) one concludes ^R(u>) = and using formula 
3.466 of Gradshteyn & Ryzhik (2000) one obtains from its 
real part the dispersion relation 
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where erfc denotes the complementary error function. As 
is well known perturbations will grow on length scales 
larger than the Jeans length y/n<r 2 /Gpb- However, dark 
halos are thought to be dynamically hot systems and 
their Jeans lengths will be of the order of the size of the 
halos themselves. Thus this part of the solution of the 
Boltzmann equation (3) is uninteresting in the present 
context and will be not considered in the following. 



3.2. Response of the halo to a density wave in the disk 

I concentrate now on the remaining part of the Boltzmann 
equation (3), 
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= 0, 



dx du dy dv dz dw 

which describes the halo response to a perturbation in the 
disk. If the gravitational potential perturbation of the disk 
is Fourier expanded the Fourier terms have the form (cf. 
equation 33 of Fuchs 2001) 



i(k x x+k y y)-k^\z\ 
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(17) 



k 2 + k 2 . This can be converted to the 



form as in equation (9) by introducing Fourier coefficients 

1 roc 1 £ 

Note that the coordinate y, which is defined in the refer- 
ence system of the disk, is related to the y coordinate in 
the reference system of the halo due to the motion of the 
center of the shearing sheet as 



(12) y^y-r flot. 



(19) 



Fourier transforming the distribution function of the halo 
particles as /hi = / d 3 kfhk exp i(k, x) and changing again 
1° £j C coordinates and the corresponding velocity com- 
ponents the Boltzmann equation (16) takes the form 



hk 
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(20) 



where k = ^Jk 2 ^ + k 2 . Equation (20) has been integrated 

with respect to the velocity components perpendicular to 
the £-axis. Assuming again a Gaussian velocity distribu- 
tion for the halo particles it can be solved straightforward 
as 
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The corresponding Fourier term of the density perturba- 
tion of the halo particles is given by 
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Pb 
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dv- 



(22) 



where the time-dependent term has been split off and is 
written no longer explicitely. Next the gravitational poten- 
tial associated with this density distribution is calculated 
from the Poisson equation, 



fc $hk = 47rGp hk , 



resulting in 
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(24) 
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The integral on the rhs of equation (24) is not trivial and 
has to be evaluated as a Cauchy principal value and the 
pole contribution at v — ' 



-— 4- 

k + 



where the halo 



particles travel at velocities in resonance with the Doppler- 
shiftcd phase velocity of the wave. Following the method of 
Kegel (1998) I obtain for the principal value of the integral 
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where the error function with imaginary argument is re- 
lated to Dawson's integral. Since the gravitational forces 
in equation (1) have to be taken at the midplane z = 0, 
it is necessary to convert $hk from the representation in 
k space to a mixed representation in (k|| , z) space leading 
to 
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The pole contribution to the integral in equation (24) has 
to be calculated according to Landau's rule resulting in 
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and in the mixed representation 
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The final result can be formally written as 

*hk ll (« = 0) = T(w-fc y rofio,*!||)*dk ll , (29) 

where the real and imaginary parts of T are defined by 
equations (26) and (28), respectively. Thus for any given 
frequency there is a contribution both from the non- 
resonant and the resonant halo particles. I find that in 
dimensionless form 



(w - k y r tt f 
2(k 2 +k 2 )a 2 



(30) 
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where a<i denotes the velocity dispersion of the disk stars 
and fc cr it = 27r/A cr it, the critical wave number of the disk. 
The function T is illustrated in Fig. 1 for parameters typi- 
cal for the solar neighbourhood in the Milky Way. The ve- 
locity dispersion of the halo particles has been estimated 
as <7h = roOo/\/2 like in an isothermal sphere. 
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Fig. 1. The function T(ui — k y r Qo, fcn)- The real part is 
shown in the upper panel, while the imaginary part is 
shown in the lower panel. The Doppler shifted frequency 
is given in terms of the epicyclic frequency and the wave 
number fcu in units of the critical wave number. Contours 
are given at levels of 10~ 2 Gpb/n 2 . Negative values are 
indicated as dotted lines. A Q parameter value of Q = 1.4 
and a ratio of disk to halo velocity dispersions of era : <7h = 
1 : 5 are assumed. 



4. Disk dynamics 

The halo response (29) to the perturbation in the disk has 
to be inserted into equation (1). Solving the Boltzmann 
equation (1) is greatly facilitated by the fact that its form 
is identical to the case of an isolated shearing sheet with 
the replacement 



$dk - (1 + T)$ dk 



(31) 
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and I can use directly the results of Fuchs (2001) even 
though the Boltzmann equation is treated there using ac- 
tion and angle variables instead of the Cartesian coordi- 
nates as in equation (1). In particular the factor 1 + T 
is carried straightforward through to the fundamental 
Volterra integral equation (equation 68 of Fuchs 2001) 



/ dk x JC (fc x , k' x ) (1 + T(fc x , k' y , w))*k.,k^ (32) 

J — oo 

where the kernel K. is given by equation (67) of Fuchs 
(2001). rw describes an inhomogeneity of equation (32) 
related to an initial non-equilibrium state of the shear- 
ing sheet. Equation (32) is separating in the circumferen- 
tial wave number k' y . In equation (32) the wave numbers 
are expressed in units of the critical wave number fc cr it 
(cf. equation 30). This implies that the Volterra equation 
describing a shearing sheet embedded in a rigid halo po- 
tential is formally the same as that of an isolated shear- 
ing sheet, because in this case T = and the halo mass 
affects only the numerical values of the critical wave num- 
ber fc cr it and the stability parameter Q. It is advanta- 
geous to consider equation (32) transformed back from 
frequency to time domain. Splitting off the w-dependent 
term exp iui ^^h* from the kernel and making use of the 
convolution theorem of the Fourier transform of products 
of two functions leads to 



(33) 
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where the operator T denotes the Fourier transform from 
lo to time domain. In equation (33) I have assumed an 
initial perturbation of the disk at time t = so that 
$k x ,k^,t'<o = 0. The Fourier transform T can be calcu- 
lated analytically using formulae 6.317 1 and 3.952, and 
the integrals with respect to k z using formula 3.466 of 
Gradshteyn & Ryzhik (2000) leading to 
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a ~§b Io° duo sin (bLo)erl{iau))e~ a w with formula 6.317. 
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t-t' + h : ' 
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in the curly bracket repre- 



sent the contribution from the non-resonant halo particles, 
while the remaining term gives the contribution by the res- 
onant halo particles. If this is inserted into equation (33) 
it takes the form 



k',t 
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Equation (35) can be integrated numerically with very 
modest numerical effort. In Fig. 2 I illustrate the response 
of the shearing sheet now embedded in a live halo to 
an initial sinusoidal perturbation of unit amplitude. For 
this purpose I use the inhomogeneity term of the Volterra 
equation 



«(fc x ,fc x )/ kx . k , in (0) 



(36) 



derived in Fuchs (2001) with / kxik /m(0) oc 5(k x -k™). The 
response of the shearing sheet to this initial impulse is 
a swing amplification event. The radial wave number fc x 
evolves as 



k x 



k x ~\~ 2A.k y t , 



(37) 



while the circumferential wave number k' y m is constant, 
which means that the wave crests swing around from lead- 
ing to trailing orientation during the amplification phase. 
Around t — 6 the amplitudes become negative which in- 
dicates that the swing amplified density wave is also os- 
cillating. As can be seen from Fig. 2 comparing the evo- 
lution of shearing sheets embedded either in a rigid halo 
potential or a live halo this characteristic behaviour of the 
density wave is not changed by the responsive halo, but 
the maximum growth factor of the amplitude of the wave 
is enhanced by a surprisingly large amount. In Fig. 2 the 
evolution of a shearing sheet embedded in a static halo 
potential has been calculated from equation (35) setting 
T = 0. The maximum growth factor is identical to within 
1% accuracy with the maximum growth factor calculated 
by Fuchs (2001) in an alternative way. In order to test 
the accuracy of the numerical solution of the full equation 
(35) I have increased the velocity dispersion of the halo 
particles to : = 1 : 13. The halo is then very stiff 
and the enhancement of the maximum growth factor of 
the amplitude of the wave by the responsive halo is only 
5.3% so that the solution of equation (35) is very similar 
to that for a density wave in a shearing sheet in a static 
halo. If I use this, which is correctly calculated, in the 
second term of equation (35) I confirm by this iterated so- 
lution the maximum growth factor of the self-consistent 
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solution of equation (35) to an accuracy better than 1%. 
The enhanced maximum growth factor of swing amplified 
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Fig. 2. Swing amplified density wave in the shearing sheet. 
The upper diagram shows the evolution in a shearing 
sheet embedded in a static halo potential triggered by 
an impulse with unit amplitude and wave vector k m = 
(—2, 0.5)fc cr ;t. Time is given in units of 0.045 epicycle pe- 
riods. The middle diagram shows the evolution of a shear- 
ing sheet embedded in a live dark halo triggered by the 
same impulse. The lower diagram shows the difference. 
The model parameters are Aj 'f2g = 0.5, Q = 1.4, a& : Uh = 
1 : 5, Gph/n 2 = 0.01, and r Q : a d = 220 : 44. 



density waves due to a responsive halo seems to be the 
equivalent of the enhanced growth of bars of stellar disks 
embedded in live dark halos seen in the numerical simula- 
tions. However, the interaction of the shearing sheet and 
the surrounding halo is not only mediated by the resonant 
halo particles, but the non-resonant halo particles play an 
important role as well. This can be demonstrated, for in- 



stance, by keeping in equation (34) only that part of the T 
function related to the resonant halo particles and solving 
equation (35) approximately by iterating the solution for 
the isolated sheet, 
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The amplification of density waves depends critically 
on the Toomre stability parameter Q. This is illustrated 
in Fig. 3 where the response of the shearing sheet to the 
same initial impulse as in the previous example is shown, 
but assuming a stability parameter of Q = 2. As can be 
seen in Fig. 3 there is neither effective amplification of 
density waves in a shearing sheet in a rigid halo potential 
or in a shearing sheet embedded in a live dark halo. 

In Fig. 4 I illustrate the maximum growth factors as 
function of the circumferential wave number fc y for shear- 
ing sheets embedded in a rigid halo potential or a live 
dark halo, respectively. As can be seen from Fig. 4 in 
both cases maximum amplification occurs for wave num- 
bers around k y = 0.5 or in the notation of Toomre (1981) 
X = fcy 1 = 2. Thus the preferred circumferential wave 
length of the density waves is 2A cr ;t. 

Following Athanassoula et al. (1987) arguments of 
swing amplification theory have been used to constrain 
the mass-to-light ratios of galactic disks. The expected 
number of spiral arms is determined by how often the az- 
imuthal wave length 27r/fc y fits onto the circular annulus, 



m = 



27TTO 
2JL 



(39) 



As shown in Fig. 4 density waves grow preferentially with 
an azimuthal wave number around fc y w 0.5fc cr ;t. This as- 
sumes a disk with a flat rotation curve, A/fl = 0.5, but 
the preferred wave length varies with the slope of the ro- 
tation curve measured by Oort's constant A (cf. Fuchs 
2001). The observed number of spiral arms of a galaxy 
taken together with the rotation curve allows then to de- 
termine the surface density of the disk. Since the wave 
number of maximum growth of the density wave ampli- 
tudes is the same for shearing sheets embedded either in a 
rigid halo potential or in a live dark halo, the presence of 
a responsive dark halo docs not influence the constraints 
on the disk mass obtained this way. 

Another constraint on the disk mass depends on the 
implied value of the Toomre stability parameter Q (Fuchs 
1999, 2003). Q must be larger than Q = 1 in order to avoid 
Jeans like instabilities in the disks (Toomre 1964). A disk 
violating this condition will evolve rapidly by fierce dy- 
namical instabilities above the threshold of Q = 1 (Fuchs 
& von Linden 1998). On the other hand, I have demon- 
strated above that the Q value must be less than Q = 2 in 
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P^JBpPp^ Fig. 4. Maximum growth factors of swing amplified den- 

rv^^-^J^^ sity waves in a shearing sheet embedded in a rigid halo 

potential (squares) and in a shearing sheet embedded in a 
live dark halo (triangles). The same model parameters as 
above have been adopted. 
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Fig. 3. Same as in Fig. 2, but adopting Q = 2. 



order that the disk can develop appreciable spiral struc- 
ture. If the radial velocity dispersions of the stars are 
known, the allowed range of Q sets constraints on the 
disk mass, which are independent from the previous ones 
(cf. equations 6 and 7). As shown in Fig. 3 this condition 
is not changed by the presence of a live dark halo. 



5. Momentum transfer to the halo 

The analogue of the angular momentum transfer from bars 
to dark halos is in the shearing sheet the transfer of linear 
momentum to the surrounding halo particles. In order to 
describe this I adapt the discussion of Dekker (1976) to the 
shearing sheet model. The acceleration of a halo particle 
is given by 



V$, 



(40) 



where I assume a gravitational potential of the form 

$(x,t) = 3?{$(x)e^ 4 }, (41) 

with complex frequency w. Thus the overall acceleration 
of the halo is given by 



< v >=- J d 3 x J d 3 vf h V$. 



(42) 



If both the gravitational potential and the distribution 
function of the halo particles are Fourier expanded 2 , 



< v > 



- / d 3 x 



d 3 v J d 3 kik^(t)e 1 ^ 
J d 3 fc'/hk'(t)e l(k '' x) (43) 
(2ir) 3 Jd 3 v J d 3 fc*k$_ k (t)A k (i) . 



In order that the gravitational potential $ be a real 
quantity its Fourier coefficients must obey the relation 
$-k(*) = an d i0r a potential of the form (41) one 

can write 



*k(t) - + $*_ k e-^**} . 



(44) 



2 It can be shown (Binney & Tremaine 1987 and refer- 
ences therein) that if the Fourier transformed potential of a 
point mass is inserted into equation (43) this leads exactly to 
Chandrasekhar's dynamical friction formula, although with a 
Coulomb logarithm defined in a slightly different way. 
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B. Fuchs: Density waves in the shearing sheet 



Again spatial coordinates £, 77, and £ are introduced with 
the £-axis parallel to the wave vector k and using equation 
(11) one obtains from equation (43) 



< v > 



i(2ir) 3 / d 3 k 



, Pb 
dv = e 
2ira 3 



} 



(45) 



x {-fc$ k k$\ — } 

1 uo + kv ^-uo*+kv > 



In equation (45) I have integrated already over the velocity 
components perpendicular to the £-axis. Two of the four 
terms of the integrand cancel out leading to 



<■£>> = -2n 3 i J d 3 k 

x fc 2 e- 29 "*|$ k | 2 { 



, Pb 

dv e 

2ira 3 



(46) 



lu + kv uj* + kv 



}• 



In the limit of — 3uj — > 



w + few o>* + fcu 
and I find 



(47) 



< v > 



-4tt 4 / d 3 k 



Pb 



2na 3 



\^\ 2 (LU-k y r n ) (48) 



x exp- 



(co - fc y r ^o) 2 
2fcV 2 



where w denotes now the frequency of the wave in the 
reference frame of the shearing sheet as in the previous 
sections. Note that according to the notation used here 
(lu — fcyTo^o) < for waves traveling along the y-axis 
in the forward direction, so that < v > > for such 
waves. From the form of the Fourier coefficients (18) it 
becomes immediately clear that only planar momentum 
is transferred to the halo, < v z >= 0, as expected from 
symmetry reasons. < v x > and < v y > depend on the 
exact form of the power spectrum |$k| 2 - Equation (47) 
shows that linear momentum transfer to the halo is medi- 
ated entirely by halo particles on orbits resonant with the 
wave. This is exactly the same as for angular momentum 
transfer from bars to halos (Weinberg 1985, Weinberg & 
Katz 2002, Athanassoula 2003, Sellwood 2003). It has also 
been known for a long time for the angular momentum 
exchange between spiral density waves and galactic disks 
(Lynden-Bell & Kalnajs 1972) or for momentum transfer 
from longitudinal plasma waves (Stix 1962). It is interest- 
ing that the exchanged linear momentum scales as Pblo\ 
in equation (48). This implies that at given mass a dy- 
namical hot halo absorbs less momentum than a cooler 
one. On the other hand, the shearing sheet would loose 
less momentum. The same effect has been described by 
Athanassoula (2003) for the angular momentum transfer 
from bars to halos. 



6. Summary and conclusions 

If a self-gravitating shearing sheet is embedded in a live 
dark halo, the halo particles respond unexpectedly strong 
to density waves in the sheet. The interaction between the 
density waves and the halo particles is mediated both by 
halo particles on orbits in resonance with the waves and 
on non-resonant orbits. If the embedded shearing sheet 
is initially perturbed by a sinusoidal wave, a swing am- 
plified density wave develops in the disk, which is of the 
same type as in an isolated sheet or a sheet in a static halo 
potential, but with an amplitude enhanced by a surpris- 
ingly large amount. This appears to be the equivalent of 
the enhanced bar growth in stellar disks embedded in live 
dark halos instead of static halo potentials seen in numer- 
ical simulations. There is transfer of linear momentum of 
density waves in the shearing sheet to the halo particles. 
This is mediated, however, entirely by halo particles on 
orbits in resonance with the waves similar to the torque 
exerted by bars on the surrounding halo. 

Acknowledgements. I thank E. Athanassoula for stimulating 
discussions and the anonymous referee for helpful comments. 
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